Identification and Optimization of RNA-Splicing Modulators as Huntingtin Protein-Lowering Agents for the Treatment of Huntington's Disease.

Huntington's disease (HD) is caused by an expanded CAG trinucleotide repeat in exon 1 of the huntingtin (HTT) gene. We report the design of a series of HTT pre-mRNA splicing modulators that lower huntingtin (HTT) protein, including the toxic mutant huntingtin (mHTT), by promoting insertion of a pseudoexon containing a premature termination codon at the exon 49-50 junction. The resulting transcript undergoes nonsense-mediated decay, leading to a reduction of HTT mRNA transcripts and protein levels. The starting benzamide core was modified to pyrazine amide and further optimized to give a potent, CNS-penetrant, and orally bioavailable HTT-splicing modulator 27. This compound reduced canonical splicing of the HTT RNA exon 49-50 and demonstrated significant HTT-lowering in both human HD stem cells and mouse BACHD models. Compound 27 is a structurally diverse HTT-splicing modulator that may help understand the mechanism of adverse effects such as peripheral neuropathy associated with branaplam.

Correction to Identification of Quinoline-Based RIP2 Kinase Inhibitors with an Improved Therapeutic Index to the hERG Ion Channel.

[This corrects the article DOI: 10.1021/acsmedchemlett.8b00344.].

Identification of Quinoline-Based RIP2 Kinase Inhibitors with an Improved Therapeutic Index to the hERG Ion Channel.

RIP2 kinase was recently identified as a therapeutic target for a variety of autoimmune diseases. We have reported previously a selective 4-aminoquinoline-based RIP2 inhibitor GSK583 and demonstrated its effectiveness in blocking downstream NOD2 signaling in cellular models, rodent in vivo models, and human ex vivo disease models. While this tool compound was valuable in validating the biological pathway, it suffered from activity at the hERG ion channel and a poor PK/PD profile thereby limiting progression of this analog. Herein, we detail our efforts to improve both this off-target liability as well as the PK/PD profile of this series of inhibitors through modulation of lipophilicity and strengthening hinge binding ability. These efforts have led to inhibitor 7, which possesses high binding affinity for the ATP pocket of RIP2 (IC50 = 1 nM) and inhibition of downstream cytokine production in human whole blood (IC50 = 10 nM) with reduced hERG activity (14 µM).

Biopharmaceutical Informatics: supporting biologic drug development via molecular modelling and informatics.

OBJECTIVES: The purpose of this article is to introduce an emerging field called 'Biopharmaceutical Informatics'. It describes how tools from Information technology and Molecular Biophysics can be adapted, developed and gainfully employed in discovery and development of biologic drugs. KEY FINDINGS: The findings described here are based on literature surveys and the authors' collective experiences in the field of biologic drug product development. A strategic framework to forecast early the hurdles faced during drug product development is weaved together and elucidated using chemical degradation as an example. Efficiency of translating biologic drug discoveries into drug products can be significantly improved by combining learnings from experimental biophysical and analytical data on the drug candidates with molecular properties computed from their sequences and structures via molecular modeling and simulations. SUMMARY: Biopharmaceutical Informatics seeks to promote applications of computational tools towards discovery and development of biologic drugs. When fully implemented, industry-wide, it will enable rapid materials-free developability assessments of biologic drug candidates at early stages as well as streamline drug product development activities such as commercial scale production, purification, formulation, analytical characterization, safety and in vivo performance.

Quantifying the Risks of Asparagine Deamidation and Aspartate Isomerization in Biopharmaceuticals by Computing Reaction Free-Energy Surfaces.

Early identification of asparagine deamidation and aspartate isomerization degradation sites can facilitate the successful development of biopharmaceuticals. Several knowledge-based models have been proposed to assess these degradation risks. In this study, we propose a physics-based approach to identify the degradation sites on the basis of the free-energy barriers along the prechemical conformational step and the chemical reaction pathway. These contributions are estimated from classical and quantum mechanics/molecular mechanics molecular dynamics simulations. The computed barriers are compared to those for reference reactions in water within GNG and GDG sequence motifs in peptides (which demonstrate the highest degradation rates). Two major factors decreasing the degradation rates relative to the reference reactions are steric hindrance toward accessing reactive conformations and replacement of water by less polar side chains in the solvation shell of transition states. Among the potential degradation sites in the complementarity-determining region of trastuzumab and between two DK sites in glial cell-derived neurotropic factor, this method identified N30T, N55G, D102G, and D95K, respectively, in agreement with experiments. This approach can be incorporated in early computational screening of chemical degradation sites in biopharmaceuticals.

Computing the Free Energy Barriers for Less by Sampling with a Coarse Reference Potential while Retaining Accuracy of the Target Fine Model.

Proposed in this contribution is a protocol for calculating fine-physics (e.g., ab initio QM/MM) free-energy surfaces at a high level of accuracy locally (e.g., only at reactants and at the transition state for computing the activation barrier) from targeted fine-physics sampling and extensive exploratory coarse-physics sampling. The full free-energy surface is still computed but at a lower level of accuracy from coarse-physics sampling. The method is analytically derived in terms of the umbrella sampling and the free-energy perturbation methods which are combined with the thermodynamic cycle and the targeted sampling strategy of the paradynamics approach. The algorithm starts by computing low-accuracy fine-physics free-energy surfaces from the coarse-physics sampling in order to identify the reaction path and to select regions for targeted sampling. Thus, the algorithm does not rely on the coarse-physics minimum free-energy reaction path. Next, segments of high-accuracy free-energy surface are computed locally at selected regions from the targeted fine-physics sampling and are positioned relative to the coarse-physics free-energy shifts. The positioning is done by averaging the free-energy perturbations computed with multistep linear response approximation method. This method is analytically shown to provide results of the thermodynamic integration and the free-energy interpolation methods, while being extremely simple in implementation. Incorporating the metadynamics sampling to the algorithm is also briefly outlined. The application is demonstrated by calculating the B3LYP//6-31G*/MM free-energy barrier for an enzymatic reaction using a semiempirical PM6/MM reference potential. These modifications allow computing the activation free energies at a significantly reduced computational cost but at the same level of accuracy compared to computing full potential of mean force.

Quantitative exploration of the molecular origin of the activation of GTPase.

GTPases play a major role in cellular processes, and gaining quantitative understanding of their activation demands reliable free energy surfaces of the relevant mechanistic paths in solution, as well as the interpolation of this information to GTPases. Recently, we generated ab initio quantum mechanical/molecular mechanical free energy surfaces for the hydrolysis of phosphate monoesters in solution, establishing quantitatively that the barrier for the reactions with a proton transfer (PT) step from a single attacking water (1 W) is higher than the one where the PT is assisted by a second water (2 W). The implication of this finding on the activation of GTPases is quantified here, by using the ab initio solution surfaces to calibrate empirical valence bond surfaces and then exploring the origin of the activation effect. It is found that, although the 2 W PT path is a new element, this step is not rate determining, and the catalytic effect is actually due to the electrostatic stabilization of the pre-PT transition state and the subsequent plateau. Thus, the electrostatic catalytic effect found in our previous studies of the Ras GTPase activating protein (RasGAP) and the elongation factor-Tu (EF-Tu) with a 1 W mechanism is still valid for the 2 W path. Furthermore, as found before, the corresponding activation appears to involve a major allosteric effect. Overall, we believe that our finding is general to both GTPases and ATPases. In addition to the biologically relevant finding, we also provide a critical discussion of the requirements from reliable surfaces for enzymatic reactions.

Quantifying the mechanism of phosphate monoester hydrolysis in aqueous solution by evaluating the relevant ab initio QM/MM free-energy surfaces.

Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).

Addressing open questions about phosphate hydrolysis pathways by careful free energy mapping.

The nature and mechanism of phosphate hydrolysis reactions are of great interest in view of the crucial role of these reactions in key biological processes. Although it is becoming clearer that the ultimate way of resolving mechanistic controversies must involve reliable theoretical studies, it is not widely realized that such studies cannot be performed at present by using most existing automated ways and that only careful systematic studies can lead to meaningful conclusions. The present work clarifies the above point by considering the hydrolysis of phosphate monoesters. The clarification starts by defining the actual issues that should be addressed in careful studies and by highlighting the problems with studies that ignore the need for unique mechanistic definitions (e.g., works that confuse associative and dissociative pathways). We then focus on the analysis of the proton transfer (PT) pathways in phosphate hydrolysis and on recent suggestions that PT involves more than one water molecule. Here we point out that most of the studies that found a proton transfer through several water molecules have not involved a sufficient systematic search of the relevant reaction coordinates. This includes both energy minimization approaches as well as a recent metadynamics (MTD) simulation study. To illustrate the crucial need of exploring the potential surfaces reliably, rather than relying on automated approaches, we present here a very careful study of the free energy landscape along a 3D reaction coordinate (RC) exploring both the standard 2D RC, comprised of the attacking and leaving group reaction coordinates, as well as of the proton transfer (PT) coordinate. Our study points out that QM/MM minimization or MTD studies that concluded that the hydrolysis of phosphate monoesters involves a PT through several water molecules, have not explored carefully the single water (1W) path (that involves a direct PT form the attacking water molecule to the phosphate oxygen). Furthermore, we identified the most likely reason for the difficulty in finding the 1W path by QM/MM minimization methods, as well as by the current MTD simulations. We also discuss the problems with current studies that challenge the phosphate as a base mechanism and emphasize that all recent studies found associative/concerted paths (although many have not realized the meaning of their results). Finally, although we clearly do not have the last word about the 1W versus 2W paths we believe that we illustrated that the crucial mechanistic problems with alternative pathways should not be resolved by just running black box search approaches.

Exploring, refining, and validating the paradynamics QM/MM sampling.

The performance of the paradynamics (PD) reference potential approach in QM/MM calculations is examined. It is also clarified that, in contrast to some possible misunderstandings, this approach provides a rigorous strategy for QM/MM free energy calculations. In particular, the PD approach provides a gradual and controlled way of improving the evaluation of the free energy perturbation associated with moving from the EVB reference potential to the target QM/MM surface. This is achieved by moving from the linear response approximation to the full free energy perturbation approach in evaluating the free energy changes. We also present a systematic way of improving the reference potential by using Gaussian-based correction potentials along a reaction coordinate. In parallel, we review other recent adaptations of the reference potential approach, emphasizing and demonstrating the advantage of using the EVB potential as a reference potential, relative to semiempirical QM/MM molecular orbital potentials. We also compare the PD results to those obtained by direct calculations of the potentials of the mean force (PMF). Additionally, we propose a way of accelerating the PMF calculations by using Gaussian-based negative potentials along the reaction coordinate (which are also used in the PD refinement). Finally, we discuss performance of the PD and the metadynamics approaches in ab initio QM/MM calculations and emphasize the advantage of using the PD approach.

Paradynamics: an effective and reliable model for ab initio QM/MM free-energy calculations and related tasks.

Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the S(N)2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.

Molecular thermodynamic modeling of the morphology transitions in a solution of a diblock copolymer containing a weak polyelectrolyte chain.

For a solution of the diblock copolymer composed of a hydrophobic block and a weak polyelectrolyte block, we obtain regions of stable aggregate morphologies in pH-solution salinity plane with the aid of the self-consistent field theory in the strong-segregation approximation. Lamellar, cylindrical, branched cylindrical, and spherical aggregates have been considered in the large interval of pH and salinity. The morphology stability maps are obtained to help control self-assembly of aggregates by variation of pH and salinity of the medium. In qualitative agreement with experiment, our calculations predict the coexistence of long wormlike micelles with branched and spherical micelles in transition zones. We compare the results of our calculations with available computer simulation and experimental data on micelles and brushes (planar and curved) formed by a diblock copolymer with one polyelectrolyte block. We show that for both weak and strong polyelectrolytes the agreement between the theory and experiment is satisfactory in most systems.